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(S| ■ Abstract 

<n: 

Q. . Within the Standard Model, and if one assumes that soft rescattering effects 

are negligible, the CP asymmetry Aq P (B — > it K) is predicted to be very 

C^; small and the ratio R = BR(B d -* tt t K ± )/BR(B ± -> ^K) provides a 

O ■ bound on the angle 7 of the unitarity triangle, sin j < R. We estimate the 

corrections from soft rescattering effects using an approach based on Regge 

phenomenology, and find effects of order 10% with large uncertainties. In 

.£h ! particular, we conclude that »4qp ~ 0.2 and sin 2 7 ~ 1.2R could not be taken 

/S ' unambiguously to signal New Physics. Using SU(3) relations, we suggest 

C$ • experimental tests that could constrain the size of the soft rescattering effects 

thus reducing the related uncertainty. Finally, we study the effect of various 

models of New Physics on A^p and on R. 
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I. INTRODUCTION 

Heavy quark decays serve as a powerful tool for testing the Standard Model and provide 
invaluable possibilities to study CP violation. However, the interpretation of experimental 
observables in terms of fundamental parameters is often less than clear. Rare hadronic decays 
of B mesons, for example, proceed through both tree level Cabibbo-suppressed amplitudes 
and through one loop penguin amplitudes. On the one hand, this situation allows direct 
CP violating effects that may give the first evidence for CP violation outside the neutral 
kaon system. On the other, these competing contributions complicate the extraction of 
Cabibbo-Kobayashi-Maskawa (CKM) angles and, in particular, the angle 7, 



7 = arj; 
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The CLEO collaboration has presented combined branching ratios for B^ — > tt ± K and 
Bd — > n Zf K ± JlJ making these modes of particular interest. In the Standard Model, these 
decays are mediated by the AB = 1 Hamiltonian, which takes the form 
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+ H.c. . (1.2) 



The flavor structures of the current-current, QCD penguin, and electroweak penguin op- 
erators are, respectively, Qf 2 ~ sqqb, Q s z 6 ~ sb^tfq', and Qj 10 ~ sbJ2e q rq'q', where 
the sum is over light quark flavors 0. The Wilson coefficients Cj are renormalization scale 
dependent; at a low scale /i ~ 1 GeV, they roughly satisfy Ci j2 = 0(1), C?,^^ = O(10~ 2 ), 
and 6*7^,10 < O(10~ 3 ). In many extensions of the Standard Model, the effect of New Physics 
at low energies is simply to modify the values of the Wilson coefficients. 

A scale-independent way of decomposing the decay amplitudes of interest is to do so 
not naively according to contributions of the operators Qf , but rather according to their 
dependence on the elements of the CKM matrix, 

A(B+ -> vr+K ) = A+ - A+ s e«e«+, A(B_ - ^W) = A+ - A+ S e^ e i5 + , 
A(B° - n-K + ) = Al - A° us ^e iSo , A(B° -> n + R-) = A° cs - A° us e^e iSo , (1.3) 

where Sq and d~ + are CP-conserving phases induced by the strong interaction, and the depen- 
dence on the CKM phases is shown explicitly. The first and second terms in each amplitude 
correspond to matrix elements of the first and second terms in TC e ff (or their Hermitian 
conjugates), respectively. Note that each term is by itself scheme and renormalization scale 
independent. 

We will avoid, as much as possible, the common terminology of "tree" versus "penguin" 
contributions, which can lead to much unnecessary confusion. The standard convention is 
to take "tree" contributions to a given decay to be those which are mediated by the current- 
current operators Q\ s 2 , and "QCD penguin" contributions to be those mediated by Qs S __q- 
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However, this is not a scale-invariant decomposition. If the computation of physical matrix 
elements could be accomplished perturbatively, this would not be such a serious failing, 
as the scale-dependence of the matrix elements would cancel explicitly against that of the 
Wilson coefficients. But this is not the case; rather, the matrix elements must be modeled 
phenomenologically, and the manifest scale-independence of the result is lost. The greatest 
difficulty is found when one is interested, as we will be, in a significant contribution to 
some process which arises essentially at long distances, where the the physics is intrinsically 
nonperturbative. There one can dispute endlessly, and pointlessly, about whether what one 
is computing is "really" tree or penguin in nature. The question obviously has no unique 
answer, but its resolution is, fortunately, of no practical consequence. 

That having been said, one can still make some general statements about the expected 
relative contributions of the operators in 7i e R to a given exclusive decay mode. The elec- 
troweak penguin operators are commonly neglected, since the contributions with a sizable 
Wilson coefficient, C9Q9, are color suppressed or require rescattering from intermediate 
states. In this case isospin symmetry of the strong interactions leads to the simplifica- 
tion A° cs = A+. It is now believed that the current- current operator contributions to 
A®' s + are roughly of same order as the QCD penguin operator contributions [0,0. To 
be specific, this observation is based on the fact that at next-to-leading order in QCD 
perturbation theory, it holds for the corresponding parton model decays b — > uus and 
b — > dds. The contribution of the current-current operators to A^ s is also expected to be 
of the same order, despite the CKM suppression, because of the large value of C2, namely, 
V u bV* s C<i ~ V c bV* s C^. .6- However, since for B ± — >■ ^K the relevant quark transition is 
b — > dds, one might expect the size of A+ s relative to A+ s to be highly suppressed by the 
small ratio iKi&KTs/KfcKsl ~ 0.02. This would hold equally for the current-current and pen- 
guin operators. If, indeed, r + = A+ s /A+. ~ iKt&KTs/^c&Ksl * s a g°°d approximation, then 
there are two important consequences: 

(z) Direct CP violation could be observed, in principle, through the CP asymmetry 
A d j P =A d SUB+^n+K°), 

dir BR(B+ -> tt+K°) - BR(B- ->■ ir-~K°) _ 2r + sin 7 sin 5 + 

CP ~ BR(B+ ->• tt+K°) + BR(B~ ->• ir-W) ~ 1 - 2r + cos 7 cos 5+ + r\ ' ^ ' ' 

However, it would be small, 

A^{B + ^7i + K )<O{X 2 ), (1.5) 

where A ~ 0.22 is the Wolfenstein parameter. For example, "hard" final-state interaction 

estimates [§-§], where the u quarks in Q\\ are treated as a perturbative loop, give ^p ~ 1%. 

(ii) Mo del- independent bounds could be obtained for the angle 7 using only the combined 

branching ratios BR{B^ — > ir K) and BR(Bd — > ^K^) p|,|T0||. One can construct the ratio 



_ BR(B°^7i-K + ) + BR(B°^7i + K-) f A°\ 2 1 - 2r cos7cos5 + r 2 

K = ^ttttt: x^m , w^ttt r^r= -n: 1 — n„ s — rj) ^ 1 - b > 



BR(B+ ^7i+K°) + BR(B- ^7i-K°) \A+J 1 - 2r + cos 7 cos 5 + + r 



where r = A^ S /A CS . Assuming that A+ s and electroweak penguin operator contributions 
are negligible, the ratio ( |1.6| ) takes the simple form 



R = 1 — 2r cos 7 cos <5o + fg . (1.7) 

The observable R may be minimized with respect to the unknown hadronic parameter r , 
yielding the inequality 

R > 1 - cos 2 7 cos 2 S . (1.8) 

Since cos 2 Sq < 1, this leads to the bound 

sin 2 7< J R. (1.9) 

The bound including electroweak penguin operators is obtained by substituting 

R - R(A+ /A° cs f . (1.10) 



It has recently been observed that the modified bound could differ by as much as ±10% JIT] . 
If true, a stringent bound on 7 would be obtained if the experimental errors in the presently 
reported R exp = 0.65 ± 0.40 [|l|] were to be reduced, with the central value unchanged. 

The B — ► iiK transitions are suppressed in the Standard Model by either CKM matrix 
elements or small Wilson coefficients. As a consequence, these decays are potentially sensi- 
tive to New Physics. In particular, in the presence of New Physics, a large CP asymmetry 
can be induced, thus violating the bound (|1.5|), and R can be modified in a way that violates 



the bound (|1.9|) . The analyses leading to (|1.5|) and to ( |i~9|) , however, explicitly assume that 



the CKM angle 7 does not enter the theoretical expression for the charged decay amplitudes 
( |1.3|) . As it is usually expressed, one requires the absence of significant contributions from 
the current-current operators Q\ s 2 to this decay channel. As noted above, this assumption is 
based on the observation that the quark level decay b — > dds is not mediated direcly by Q" %■ 
However, this naive treatment of the dynamics ignores the effects of soft rescattering effects 
at long distances, which can include the exchange of global quantum numbers such as charge 
and strangeness. It is the purpose of the next section to discuss an explicit model for such a 
rescattering, and to consider its effect on the bounds fll.5| ) and (|1.9| ). Before presenting the 
model, however, we must make some general comments about what we do, and what we do 
not, expect to accomplish with this exercise. 

We will consider the contributions to B + — > tt + K° from rescattering through coupled 
channels such as B + — > n (ri)K + or corresponding multi-body decays. In fact, we will treat 
only the two body intermediate states explicitly. The model we will employ will be based 
on Regge phenomenology, including the exchange of the p trajectory and others related by 
ST/ (3) flavor symmetry. The coupling of the trajectory to the final state will be extracted 
from data on np and pp scattering cross sections. A few points are in order. First, we have 
little confidence in the quantitative predictions of the model per se. In fact, we believe that 
it would be irresponsible to claim an accuracy of better than a factor of two for the size 
of soft rescattering effects, using any model currently available. Neither our nor any other 
model should be taken as a canonical framework for the estimate of final state interactions in 
B decays. Rather, the purpose of our calculation is to be illustrative: our model will predict 
r + at the level of ten percent, with no fine tuning or unnatural enhancements. We will use 
this result to argue that such a value of r + is entirely generic within the Standard Model. 



At such a level, the effect of final state interactions on total branching fractions is small, 
but the effect on quantities which are most interesting for r + = 0, such as AqpIB^ — > ir^K) 
and R, can be more dramatic. 

Second, the rescattering effects which we will consider are not already included in an 
analysis of the strong "BSS" phases induced when the virtual particles in a penguin loop 



go onto the mass shell [JL2]. The fact that such an analysis typically yields small corrections 
cannot be used to argue that rescattering contributions to r + must be small. On the contrary, 
we will consider intermediate states with on-shell pseudoscalar mesons, rather than on-shell 
quarks. In the absence of an argument that parton-hadron duality should hold in exclusive 
processes involving pions and kaons (for which there is scant evidence), one must conclude 
that the long distance physics of meson rescattering is not probed by the BSS analysis. 
(The recent proposal that rescattering effects must be small [|13[] does not go outside the 
BSS framework in obtaining quantitative estimates.) 

Third, the issue of whether the processes we consider are of the "tree" or "penguin" type 
is a dangerous red herring. As discussed above, the question has no scale-invariant meaning, 
and also no important implications. Our model addresses contributions to the well-defined 
amplitude A£ 8 . Within the model, we will first use the current-current operators Qi% to 
generate the transition B — > 7r°(n)K~, which will then rescatter to n~K°. However, we 
will deliberately refrain from referring to this as a "tree" contribution, in order to avoid 
unnecessary confusion. 

Fourth, if one decomposes the matrix element for B~ — ► n~K into amplitudes of definite 
isospin, the rescattering process which we will consider is non-trivially embedded in their 
sum. Recent isospin analysis of this decay J13|]l4]l have stressed the importance of the 
strong phases associated with the isospin amplitudes. (For previous isospin analyses of 
B — > nK decays, see [jl5j[l6| .) In the isospin language the magnitude of A+, depends on the 
differences between these phases and vanishes in the limit that they vanish. However, the 
isospin decomposition sheds no light on the sizes of these phase differences, hence on the 
magnitude of A+ s , leaving the former as free parameters. The literature presently contains 
a variety of quite divergent opinions concerning the "natural" size of rescattering effects in 
this decay pyTT|-|il|]. The purpose of our calculation, however crude, is to address the issue 
more quantitatively by providing the first analysis to model the magnitude of rescattering in 
this decay. Finally, we note that the rescattering in question is inelastic, despite its quasi- 
elastic kinematics, and cannot be studied adequately in any model of purely elastic final 
state phases. 

The calculation of the Standard Model predictions is given in section 2. First, we describe 
how the final state interactions (FSI) affect the CP asymmetry and the bound on 7. Then 
we calculate, using the phenomenological Regge model, FSI corrections for specific two body 
states. In section 3 we suggest experimental tests that could potentially give an upper bound 
on these contributions, independent of hadronic models for the final state rescattering. In 
section 4 we analyze which types of New Physics models can significantly affect the relevant 
B — > 71K decays, and whether there are relations between such new contributions to the 
charged and neutral modes. We summarize our results in section 5. 



II. FINAL STATE INTERACTIONS 

A. The effects of FSI Corrections 

We would like to investigate the impact of final state rescattering on the CP asymmetry 
Aqp(B — > ir^K) and the ratio of branching fractions R. The rescattering process involves 
an intermediate on-shell state X, such that B — > X — > Kit. In particular, we assume 
that there exists a generic (multibody) state Kmr. In a straightforward generalization of 
Eq. ( |1.3| ), the charged and neutral channel amplitudes can be written as 

A{B + - Knrr) = A^ - ^ S W^ , (2.1) 

A(B° -> Kmx) = Af s - A^e* V 5 ° . 

Rescattering contributions, again decomposed according to their dependence on CKM fac- 
tors, are given by 
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A{B + -+ Knn -► n + K°) = S?A^+ - S^A^fe 

A(B° -> Kmx -> iT K+) = S™Af s - S^A^e^ , (2.2) 

where Sf 1 is the complex amplitude for rescattering from a given multibody final state to the 
channel of interest. Analogous contributions arise in the conjugated channels. In the limit 
where one neglects electroweak penguin operator contributions, isospin symmetry requires 
A+ s = A Q cs1 and this equality is not spoiled by rescattering effects. The i = 1,3,4 rescattering 
amplitudes can, for our purposes, be absorbed into the unknown amplitudes in Eq. (|1.3| ). 

We are interested in the possibility that the rescattering of transitions mediated by Q" % 
are significant enough to dominate A+ s , so we make the approximation 

A+ S e i5+ = E S n 2 A n u t , (2.3) 

n 

and define e = A^ s /A^ s . Let us assume that rescattering effects do not dominate the overall 
decay, so we may retain just terms linear in e. In that case, Acp of Eq. ( |1.4| ) and R of 
Eq. ( |1.6|) take the form 

dir 26 sin 7 sin £+ 

1 — 2ecos7coso + 

= 1 - 2r cos 7 cos ^ + rg 
1 — 2e cos 7 cos S + 

Once again, we may extremize R with respect to the unknown r , 

R> L^ggWj , (2 . 6) 

1 — 2ecos7cos^) + 

Following the same line of reasoning as before with respect to the unknown strong phases 
£ and S + , we find the new bound 



sin 2 7<i?(l + 2e V / l-i?), (2.7) 

or solving for cos 7, 



|cos7| >Vl-R-eR. (2.8) 

It is clear that even a small rescattering amplitude e ~ 0.1 could induce a significant shift 
in the bound on 7 deduced from R, effectively diminishing the model-independent bound of 
Ref. as R —>■ 1. It is also clear from (|2.4j ) that a small rescattering effect, again e ~ 0.1, 
could in principle generate an 0(10%) CP asymmetry which is significantly larger than the 
bound ( |1.5| ) on Aqp(B ± — * n ± K). Therefore, in order to understand whether a large CP 
asymmetry signals New Physics, and whether it is possible to obtain a bound on 7, it is 
useful to employ a particular model of soft FSI to obtain an order of magnitude estimate of 
the effect. In the next section we will estimate the amplitude e using a phenomenological 
approach based on the exchange of Regge trajectories |iT|-|T9"|| . 



B. The two body rescattering contribution 

We will estimate the contribution to e from the rescattering of certain two body inter- 



mediate states. While these channels alone are not expected to dominate rescattering |T7 
we might expect them to provide a conservative lower bound on the size of the effect. At 
any rate, it would be peculiar for the total effect of rescatttering to be significantly smaller 
than the two body contribution. The most important channels in charged B decays that 
might rescatter to the final states of interest are B~ — > tt°K~ and B~ — > r]K~. (To be 
conservative, we will neglect the rj'K channel, which is unrelated to the others by SU(3) 
symmetry. Its inclusion would likely enhance the effect which we will find.) To estimate the 
contribution to e of these channels, it is necessary to estimate both the relative amplitude 
Af^ /A£ s for producing the intermediate state, and the amplitude e<i = IS2I for it to rescatter 
to the final state of interest. In our model, then, e = 62(^+^4+), and we will extract only 
the magnitude e. We will not attempt to predict the strong phase 5 + , which is even more 
model-dependent. 



The amplitude ratio A^, /A+ s may be estimated using factorization |2~0|-|2~3| and the Bauer- 



Stech-Wirbel [[24] model, starting from 7i e //. We will ignore the electroweak penguin op- 
erators. In addition to the soft FSI contributions which we will estimate using the leading 
order Wilson coefficients, there are also "hard" FSI phases which can be generated via quark 



rescattering |12| and are estimated at next-to- leading order || . As discussed earlier, the final 
state interactions which we consider are distinct from these "BSS" phases. 

We will impose SU(3) symmetry on all aspects of our estimate of e. Corrections to the 
SU(3) limit are typically at the level of 30%, small compared with other uncertainties in 
the calculation. Hence we must consider both the intermediate states n°K + and t]K + . For 
normalization, we will also need to compute A^ s , which in the BSW model is induced by the 
QCD penguin operators. The computation is straightforward, and we find 

A 2 +(B + - P°K + ) = - G F m||K 6 lCl((Ci + C 2 /3) F^f^(m 2 p0 )L K (f, P o) 



+ (C 2 + d/3) F^ (m^L^po) | 

A+(B+ - vr+K ) = - G F m 2 B \V cb V c * s \l(C A + C 3 /3) F s K d ofl{m 2 K )L K {^) 

2m 2 1 

+ (C 6 + C 5 /3) F^o/-(m^)^M.(^) , (2.9) 

m s m 6 J 

for P° = 7i°, i]. The form factors f+(q 2 ), the decay constants Fp 192 , and the kinematic 
functions Lp(ni) and Mp(/Xj) are defined as follows: 

(Plfrf&lB) = /£Gz 2 )(pb + pp) m + / p (g 2 )fe -pp^ , 

(F 3 |gi7 M 7 5 g 2 |0) = -2v / 2Fr 2 Pp, 

/ P (g 2 ) 
Lp(^) = 1 - fi P + n , 



Mp(fi) = \ 



(3-y + (l-3y)/i P -(l-y)/i) 
/+(? ) 



(2.10) 



where /Zj = m 2 /m 2 B , q 2 is the momentum carried away by the current, and y ~ 1/2 — 1 is 
related to the distribution of quark momenta in the pseudoscalar mesons. Note that in the 
BSW model, the dominant contribution to A 2 l ~^(B + — > P°K + ) is from the current-current 
operators Qi%, while the dominant contribution to A 2 f(B + — ► tt + K°) is from the QCD 
penguins Q^ s 6 . 

The ratios (A^+/A+) 7r o r) simplify significantly if SU(3) symmetry is imposed on the 
quantities which appear in Eq. ( |2.9|) . In particular, we make use of the SU(3) relations 

FJ = Fg/y/2, ff(0) = f?(p)/y/2, 

F^ = Ff/VE 1 fl(0) = f«(0)/VQ. (2.11) 

We have checked that the inclusion of 77 — 77' mixing, which violates £77(3), would change 
our final answer by no more than 20%. The ratios of interest then may be written 

fA 2 u +\ __ \V ub V: s \ (l + l/3)(C 1 + C 2 )/V2 



A+ J^ \V*V£\ (C 4 + C 3 /3) + (C 6 + C 5 /3) (2m 2 K /m s m b )(M 7T ( f i K )/L n ( f i K )) ' 

fA 2 +\ 1 (A 2 +\ 

(2.12) 



A+ L VS \A+ 



where 77 — 77' mixing is neglected. With values for the coefficients Cj taken from Ref. p5 |, at 
leading order and with A^ = 225 MeV, we find \(Al+/A+)^o\ ~ 0.35. 

We now turn to an estimate of the rescattering amplitude e 2 . Our technique is described 
in detail in Ref. |17|| , and here we only outline the procedure. We begin by writing an 
expression for the discontinuity of the amplitude for the charged B decay, 



Disc A(B~ -> K^-) FSI = - I f d ^ 1 , d l^ 2 (2n) 4 S^(p B - Pl -p 2 ) 
v ' b 2 J (2tt) 3 2Ei (2tt) 3 2£ 2 V ; V ; 

xA(B" -► 7r°(pi)ii:-(p 2 )) A1(ir (pi)iif-(p 2 ) -► X°7r-) . (2.13) 

The rescattering matrix element A / t(7r°(pi)iC — (P2) "^ -^°vr~) has a well known parameteri- 
zation inspired by Regge phenomenology [fT7| . For the exchange of the leading p trajectory, 
it may be written as 

p-Mra(t)/2 / o \ a(t) 

A< + (7r ir- -> KH~) = 7 (t) - , (2.14) 

cos(7ra(t)/2) \so/ 

where 7(f) is a residue function, s and £ are the Mandelstam variables, and so = 1 GeV 2 is an 
arbitrary hadronic scale. We take a linear Regge trajectory, a(t) = «o + oc't, with «o = 0.44 
and a' = 0.94 GeV~ 2 . Since «o is approximately the same for the p, K* and u> trajectories, it 
is convenient to consider a single "octet" trajectory which carries a representation of SU{3) 
and contains all the vector mesons. Taking j(t) = 7(0) = 7 for simplicity, the discontinuity 
can be calculated, 

s \ao-l 



where 



Disc A(B~ -^ K ir~) FSI = ^(s) ( — ) A(B~ -> ^K~) , (2.15) 



\ I e -i-n-a /2 

H s ) = TH 1 J^\ m f 1 \ ~1^\ ■ ( 2 - 16 ) 

lo7r cos(7ra;o/2j s a' (m(s/so) — m j 2) 

We restore the FSI contribution to A{B~ — > K 7r~) by use of a dispersion relation. The 
dispersion integral may be evaluated in closed form with the approximations «o = \ and 
ln(s/s ) = ln(m|/s ), 

1 r°° ds ( s\ a °~ 1 



A(B~ - KH~) FSI = 7 e 2 (m|) A(£T - Tr ^") x - / 



7T Jjm.+my) 2 S — m| \Sq 



= iye 2 (m B ) ^A(B~ - vr 1 ^) , (2.17) 

where we have taken the limit Tn^/ K \/m B ~~ * 0- For rescattering through the r]K~ channel, 
the only difference is an 577(3) group theory factor in the residue function 7. The magnitude 
of the contribution of a given channel to the soft rescattering amplitudes defined in the 
previous section is then 

e 2 = 7 v^|e 2 ( m |)|, (2.18) 

m B 

where the residue 7 depends on the channel. 

Finally, we must make a numerical estimate of 7, which parameterizes the coupling of the 
p trajectory to the pseudoscalar mesons. These couplings may be estimated by considering 
pp and 7i + p scattering data and using 577(3) relations. The exchange of the p trajectory 
requires a coupling at each vertex, so 
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lp{VV -» PP) = 7; 



pppl 



J p (np -► Tip) = Jxnplppp ■ 



(2.19) 



The optical theorem gives 



&tot 



-M 



f-*f 



^2.20) 



for both pp and Tip scattering. The forward scattering amplitude is obtained from Eq. ( (2.14|) , 
written for pp and np scattering, by setting t — (that is, cos^ = 1). The residue func- 



tions ( 2.19 ) entering the expression for the forward scattering amplitudes are fixed from the 
Particle Data Group parameterizations of pp and Tip scattering data [[2(J , 



a 



ik 
tot 



Xjk I — 

so 



0.08 



+ Y lk [~ 

SQ 



-0.56 



(2.21) 



The first term represents the Pomeron contribution. The second comes from the p trajec- 
tory for 7T + p scattering, and is a combined contribution from the p and a trajectories for pp 
scattering. The Particle Data Group fit gives Y pp = 56.08 mb and Y n + p = 27.56 mb. Assum- 
ing that the p and a trajectories contribute equally in the pp channel, and using again the 
approximation ao = |, we find 



2 _ 2 

/0 ITTTTp 



2-So^rp 



Y 



72. 



^2.22) 



/v 



As defined, 70 is the 7r + 7r~p° coupling, since of the vector meson octet, only the p contributes 
in the ir + p channel. The residue functions which we will require may be found by applying 
SU(3) symmetry. For p exchange in the nK channel, we find 7^^ = 70, ^kk p = —lo/V^, 
and T^p = 0. As a consequence, we have for the p trajectory, 



-f p (n°K- -► n-K°) 



V2 



7o 



1 p (r]K- -» n-K°) = 



and similarly for K* exchange, 







(2.23) 



(2.24) 



With these residues, we complete the computation of the rescattering amplitude €2- 

We are now in a position to estimate the contribution to e from a given two body rescat- 
tering channel, by combining the rescattering amplitudes ( 2.18 ), the production amplitude 
ratios fl2.12| ), and the residues ( |2.23| ) and ( |2.24| ). For example, for the intermediate ir°K~ 



state, rescattering via p exchange, we find a term 



-nKp 



An+ 

A+ 



m B 



e 2 (m 2 B )\ lp (7r°K~^n~K0) 



(2.25) 



We do not know the strong phase which multiplies e n K P - In adding the three contributions 
t-KKp, tnKK*, and e V KK* (recall that e V K P vanishes), we must make some assumption about 
their relative phases.^ For the purpose of this estimate we might imagine adding them 
incoherently, so e 2 = e 2 KKp + & KK * + e^Kp- Our estimates for the various channels are 
e^Kp — 0.044, e-xKK* — 0.022, and € v kk* — 0.022, which by this prescription would give 
e ~ 0.06. Alternatively, adding them coherently would yield e ~ 0.09. These unknown 
relative strong phases are one important source of uncertainty in this estimate of e. 

Of course, there are many other sources of uncertainty as well. Perhaps the most severe of 
these is the neglect of multibody intermediate states, which we have omitted because we have 
no good model for them. For this reason, we are likely to have, if anything, underestimated 
the effect of FSI. More model-dependence arises from the use of factorization to estimate 
the ratios Af^ /A+ s , given that this ansatz is known at times to fail in B decays to light 
pseudoscalar mesons. Smaller uncertainties arise from the phenomenological extraction of 
the residues and from the use of SU(3) symmetry. Finally, we note that in some (but not 
all) models of Regge exchange, there may be cancelations which further suppress the small 
contribution from the K* exchange to both the 7]K~ and tc°K~ intermediate states. (For 
recent discussions of these questions, see Refs. p7|-p9[].) Similar effects may enhance the 
contribution from 77'i^ _ .[] Our neglect of the rj' has been conservative, in the sense that 
it has likely caused us to underestimate the total effect of rescattering. We have done so 
because, in the presence of the anomaly, the rj' is not related by unitary symmetry to any 
other meson. Hence we have little guidance, other than from the quark model, for how to 
include it. In the end, it certainly would be unwise to trust our estimate of e to better 
than a factor of two, and even that much confidence would be optimistic. The same caveat 
should be applied to any phenomenological model of soft final state interactions. What is 
important here is that we have found neither a dominant effect of order one or larger, nor 
an insignificant effect of order one percent. 

C. The effect of FSI on A$p and R 

We close by returning to the effect of FSI on the observables Aq p and R. We have seen 
from our model that rescattering effects as large as e ~ 0.1 — 0.2 easily could be consistent 
with the Standard Model. Therefore, unless 7 were known independently to be very small, an 
observation of A^p ~ 0.2 could not be taken unambiguously to be a signal for New Physics. 
Note that we do not claim to predict such a large asymmetry; we simply observe that it 
would be neither unnatural nor surprising for it to be generated by final state interactions. 



*We prefer not to impose constraints from SU(3) symmetry on the phases because of the substan- 
tial model-dependence already present in our calculation. Doing so would not change substantially 
the magnitude of the rescattering effects, but it would imply a smaller uncertainty then we would 
advocate. 

'We are grateful to H. Lipkin for discussions of this point, and for stressing to us the important 
role played by tensor meson exchange. 
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On the other hand, we would maintain that a larger asymmetry, such as Aqp ~ C(l)> would 
still be an exciting sign of a source of CP violation beyond the CKM matrix. 

Similarly, we can consider the effect of e ~ 0.1 on the bounds ( |2.5| ) and ( |2.7p on 7. For 
example, the fractional correction to the bound on | 00S7I is A = eR/^1 — R. The value of 
A is a strong function of the experimentally observed i? ex p, 

A ~ e for R exp = 0.65 , , * 

A~2e for R cxp = 0.80 . l ] 

The bound deteriorates quickly as R exp — > 1. In terms of sin 2 7, we certainly would conclude 
that the observation sin 2 7 ~ 1.2R exp would not constitute an unambiguous signal of New 
Physics. 



III. MODEL INDEPENDENT BOUNDS ON THE FSI CORRECTIONS 

Our phenomenological model suggests that soft FSI contributions of the current-current 
operators Q^ s 2 could account for 0(10%) of the B + —>■ K°tt + amplitude. This has a dramatic 
consequence, namely making Acp ~ 10 — 20% realistic within the Standard Model. In 
view of the large theoretical uncertainties involved, it would be extremely useful to find 
an experimental method by which to bound the magnitude of the FSI contribution. The 
observation of a larger asymmetry would then be a signal for New Physics. In this Section 
we describe such an attempt, along the lines proposed in Ref. Q] for the decay B — >■ cf)K. 



The idea is to find decay modes mediated by the quark level transition b — > ssd, for which 
branching ratio measurements or upper bounds would, by application of flavor SU(3) flavor 
symmetry, imply a direct upper bound on e. 

The most interesting modes in our case turn out to be B ± — > K ± K. The effective 
Hamiltonian for b — > d transitions may be obtained from ( |1.2|) by the substitution of s — - > d 
in the operators and in the indices of the CKM matrix elements. In analogy with Eq. ( |1.3| ) 
the amplitudes may be decomposed according to their dependence on CKM factors, giving 

A(B + -> K + W) = A cd - A ud e l ^e lS , 

A(B~ -► K-K°) = A cd - A ud e-^e iS . (3.1) 

Invariance under the SU(3) rotation exp(i^X 7 ), i.e., interchange of s and d quark fields, 
implies equalities among operator matrix elements, 

(K-K \Qf\B~) = (K n-\QT\B-), q = u,c,; z = l,2 (3.2) 

{K-K°\Q}\B-) = (K\-\Q s t \B-), 1 = 3, .., 10 . 

These lead to the relations 

A ud e i5 = A+ s e l5 +^(l + R ud ), A cd = A+^(l + R cd ) , (3.3) 

Vus ' cs 

where R ud and R cd parameterize SU(3) violation, which is typically of the order of 20 — 30%. 
Note that it is only an SU(2) subgroup of 577(3), namely ?7-spin, which is required to derive 
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these relations. Since the B~ carries U — and the transition operators Q\ and Qf carry 
U — |, it is only the U — h component of the K~K° final state which couples to the decay 
channel. As a result, we have the freedom to add any pure U = § combination to the final 



state, involving additional tt~tt° and ir~r] pairs, without affecting the relations (|3.2|) . As it 
turns out, it is the simple combination K~K° which yields the most phenomenologically 
interesting bound. 

An upper bound on e follows from the ratio 



R 



K 



BR(B+ -► K + K°) + BR{B- -> K~K 



-vro\ 



BR(B~ 
After some algebra, we obtain 



KH+) + BR(B~ -» K 7T-) 



e < XJR K (1 + Re[i^d]) + \ 2 (Rk + 1) cos 7 cos ,5+ + C(A 3 , X 2 R ud , cd ) 



(3.4) 



(3.5) 



where A ~ 0.22 is the Wolfenstein parameter. (In deriving this relation, we assume that 
the full matrix element for B~ — > K~K° is not smaller, in magnitude, than the partial 
contribution from rescattering. It is possible, if there is some fine tuned cancelation, for 
rescattering actually to lower the branching ratio; in this case, the relation ( |3.5|) would not 
be valid.) This bound becomes more reliable if we set cos 7 cos 5 + = 1. With Eq. (|2.7|), we 
find a bound on 7, 



sin 2 7 < R + 2XRJR K (1 - R) + 0(X\ XR ud , cd ) 



(3.6) 



ignoring electroweak penguin operators. One could include them, as before, with the sub- 
stitution in fll.lOp . 

We are also interested in obtaining an upper bound on A^p. Keeping cos^ free gives 



\A$p\ < 2esin 7 + e>(e d 



(3.7) 



In the absence of a nontrivial bound on 7 from Eq. (|3.6| ), Aq P is maximized at sin 7 = 1, 
leading to 

|^|,| < 2A V /^(1 + Re[R ud }) + 0(X 3 , X 2 R udfid ) . (3.8) 

However, an independent bound on 7 would place a tighter limit on the asymmetry, 



\A A &,\ < 2A y/R K R (1 + Re[R ud \) + 2A 2 v^R {R K Vl -R + R K + l) + 0{X 3 , X 2 R ud , cd ) . (3.9) 
Again, electroweak penguin operators may be included with the substitution ( |1.10 ). Follow- 



ing Ref. J30l , analogous bounds on Aq P (B ± 
Rk the ratio 



0if(*)±) can be obtained by substituting for 



/JBR(B ± -► <f>K d 



'BR{B ± -> K <.*)K±) + ^BR{B ± -> ^{p*)) 

The CLEO Collaboration has recently obtained the upper bound Rk < 0.95 at 90% 
C.L., including only statistical errors, and approximately Rk < 1-9 once systematic errors 
are included as well ||31|| . Unfortunately, this is not very restrictive, implying e < 0.4 and 
Aqp < 0.6. It is possible that more interesting constraints on e and Aq P could be obtained 



12 



in the future, given that Rk ~ 0(X 2 ) in the limit of vanishing FSI. Ultimately, the utility of 
such a bound is limited by the fact that, for e small enough, Rk becomes independent 
of e, since rescattering channels are then negligible compared with other contributions. 
Conversely, a large observed value for Rk would not necessarily mean that FSI contributions 
are large, since this could also result from New Physics which enhances the b —>■ d transitions. 
Resolving the question of how to distinguish these two possibilities is left for the future. 

IV. NEW PHYSICS 

Assuming that a large CP asymmetry is observed in B ± — > ir^K decays or that sin 7 
is measured independently and found to violate the bound Q1.9Q , this could be explained 
either by large soft rescattering effects or by New Physics. From our calculations above, 
it is clear that theory cannot at present exclude the possibility of large FSI effects. Yet, 
it is also possible that these effects are not large and that, in the future, the experimental 
tests proposed in the previous section will imply that New Physics indeed is required. It is 
important, then, to understand which extensions of the Standard Model could contribute 
significantly (and with new CP violating phases) to the relevant B — > -kK modes. Further- 
more, we would like to understand whether New Physics contributions to the charged and 
neutral modes should be expected to have any special features (such as isospin symmetry 
relations) which will allow further tests. 

The most general New Physics effects on the amplitudes (|1.3| ) can be parameterized by 
six new parameters, 

A(B° -> n~K + ) = A° SM - A^e^e^, A(W -+ 7r + K~)= A° SM - A u N e~ i(t>u e iSu , 

A(B + -> n + K°) = A^ M - A d N e i<t)d e iSd , A(B~ -► 7T#°)= A SM - A d N e~^ d e iSd . (4.1) 

Here Asm are the Standard Model amplitudes, where A and A carry the same strong 
phases and opposite weak phases, as do A ± . We introduce the New Physics amplitudes 
A u ^ , with CP violating phases (j) Ut d and strong phases 5 u ,d- 

A first class of models are those which give potentially large tree-level contributions to 
AB = 1 four quark operators. These include supersymmetry without R parity and models 
with extraquarks in vector-like representations of the SM gauge group. A second class of 
models gives new contributions to four quark operators through loopdiagrams only, with 
new particles running in the loop. This includes various supersymmetric flavor models and 
a sequential fourth generation. A third class of models gives new contributions tothe AB = 1 
chromomagnetic dipole operators. Finally, there are models where no significant effect is 
expected. These include models of extra scalars and Left-Right Symmetric models. 

A. Supersymmetry without R Parity 

In supersymmetric models without R p there are new, slepton mediated tree diagrams, 
contributing at tree level to the b — > uus and b — > dds transitions. 

First, consider the X[j k LiQjdk couplings. For b — > uus transitions, the contributions come 
from charged slepton mediated diagrams, while for b — > dds transitions, the contributions 
come from sneutrino mediated diagrams, 
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Aire** oc Y ^B 2 , (4.2) 

Ad i<f> d V^ A il2, A il2 ~r A t23 A ill T A i21 A i31 T A ill A i32 / x o\ 

* fe m*Pi) [ ' } 

We learn that indeed the new physics introduces six new independent parameters. In the 
special case where (i) m 2 (£^) = rr?{y^) and {%%) A^ 13 A^* 2 is much larger than the other three 
combinations that appear in ( |4.2| ), isospin is a good symmetry (similar to the SM QCD 
penguin diagrams). The first condition is fulfilled in many models, but the second is not. 
Generically, we have A d N ^ A%, (fid ^ 4>u an d S u ^ 5 a- 

Second, consider the \'-j k Uidjdk couplings and note that these couplings are antisymmetric 
in (j, k). For b — > uus transitions, the contributions come from the down squark mediated 
diagrams only, while for b — > dds transitions, the contributions come from up, charm and 
top squark mediated diagrams: 

Ale* oc %^f , (4.4) 

Ajre** oc V %^ . (4.5) 

Again, there are six new independent parameters unless (i) the z = 1 contribution dominates 
/L^r and {%%) m 2 (di) = m 2 (ui), in which case isospin is a good symmetry. The first condition 
is unlikely to be fulfilled. We expect A d N ^ A^, (fid ^ (fi u and 5 U ^ Sd- 

B. Singlet Down Quarks 

Models with additional S , f/(2)i-singlet down quarks could lead to dramatic effects in CP 



violation in the interference between decays with and without mixing [32]. This is a result 



of Z-mediated tree level contributions to B — B mixing. Could there also be large effects in 

yldir? 

Obviously, there are new contributions to the decay amplitudes from Z-mediated tree 
level diagrams, 

A u iK k UsbQZv* (4 6) 

m z 

A di^ d(X U^dd (47) 

™>z 

where gzff oc (T^ —Qj sin 2 6w), and U s b is the mixing angle in the Z couplings. The present 
experimental bound on BR(B — > Xfi + fi~) gives \U s b/(V c bV cs )\ < 0.04. This suppression 
roughly compensates for the loop suppression of the standard model QCD penguin. Thus, 
we learn that 

(i) A N could be of the same order as A SM ; 

(ii) A^j ^ A d N because of the gz qq factor; 

(in) (fi u = (fi d = axg(U sb ). 
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C. Supersymmetric Flavor Models 

Supersymmetric models (with R p conserved) contribute to b — > s transitions through 
penguin diagrams with squark - gluino loops. The ratio between this contribution and the 



standard model penguin amplitude may be estimated to be [33 



A SM ~ a 2 max(m?,M 3 2 ) V tb V t * s \n(m 2 t /m? c ) ' l ' ' 

where K is the mixing matrix for the gluino couplings and rj ~ (m| — r??J)/(m~ + m~) is a 
measure of the non-universality in the squark masses. We see that even in the absence of a 
super-GIM mechanism, namely with 77 ~ 1, the supersymmetric contribution is comparable 
to the Standard Model one only if the relevant mixing angle, namely K^ 2 K 33 , is large. Such 
a situation is phenomenologically allowed. Furthermore, the fact that m s /mb ~ \V c b\ implies 
that a large 5r — Sr mixing (K32 ~ C(l)) is not unlikely [31]. Below we examine whether 



this possibility is indeed realized in various supersymmetric flavor models. 

Most flavor problems of supersymmetry are solved (without giving up naturalness) in 
models where the first two squark generations are heavy [J35| -|39|l . Mixing angles with the 
third generation can be large. In these models, for bi and gluino masses in the range of 



100 - 300 GeV, it was found [^0J that the supersymmetric QCD penguins (namely, squark 
- gluino loops) can be twice as large as the standard model ones for b — > s transitions, and 
can carry a new phase. Similar to the standard model QCD penguins, isospin is a good 
symmetry here, namely A d N = A^, 4>d = 4>u and 5^ = 5 U . 

The situation is different in models of Abelian horizontal symmetries, where alignment 
of quark and squark mass matrices is the mechanism which suppresses supersymmetric con- 
tributions to flavor- changing neutral currents (FCNC) [f|l],[|^]. The constraints from K — K 



mixing require that the (32) entries in the mass matrices are also small, leading to a sup- 
pression of K 32 - A similar situation occurs in models of non- Abelian horizontal symmetries, 
where degeneracy of the first two squark generations suppresses FCNC [[36|,[43|-f|5|] . Typi- 



cally, the first two generations are in a doublet and the third in a singlet of the horizontal 
symmetry. Then, bi.R — sl,r mixing does break the horizontal symmetry and is, therefore, 
suppressed: ^32-^33 is of order V ts V t b and the supersymmetric contribution to the b —> s 
transition is typically small |33|f|6 . 



Finally, if CP is an approximate symmetry of the New Physics (which is a viable pos- 
sibility in the supersymmetric framework), so that all CP violating phases are small, say 
4>u,d — C(10~ 3 ), then we expect New Physics to contribute to A^p at the level < 10~ 3 
regardless of the size of the supersymmetric contributions to the various B — > iiK decays. 

D. Fourth Quark Generation 

Models of four quark generation require, of course, that the fourth generation neutrino 
is rather heavy (> mz/2). If this possibility is realized in nature, then we expect large 
new contributions from QCD penguin diagrams with W - t' loops. The m^-dependence of 
this contribution may compensate for a possibly small CKM factor, VfbV t * s , and become a 
significant, if not dominant, contribution. Furthermore, as the 4x4 quark mixing matrix 
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has three independent CP violating phases, this contribution is likely to carry a new phase. 
Isospin should be a good symmetry and 5 U = 5d is likely. 

E. Models with Enhanced b —> sg. 

Models with enhanced b — > sg dipole operator coefficients, leading to B(b — > sg) ~ 10%, 
have been suggested as a possible resolution of several potential puzzles in inclusive B 
decays |47| , 48| . Examples have been discussed which employ squark-gluino loops, vectorlike 



quark- neutral scalar loops, or techniscalar exchange |^SIM. Despite the large overall rate, 
the dipole induced amplitudes for rare B decays are of same order as the Standard Model 
amplitudes and so the two can interfere substantially. Arbitrary new weak phases in the 
dipole operator coefficients can therefore lead to sizable A^ P > 10% [|50| . More modest 



enhancement of 6 — > sg can also lead to sizable CP asymmetries. Such an example has been 



discussed in the context of topcolor models |5T|. Again isospin is a good symmetry. 



F. Discussion 

There exist, of course, extensions of the standard model where large new contributions 
to the relevant B — > -kK decays are unlikely. Below we give a few examples. 

Charged Higgs mediated tree diagrams contribute only to A^, so A d N = 0. The CKM 
combinations are similar to those in the standard model W mediated tree diagrams, namely 
we still have <p u = 7. The contribution is, however, small compared to the Standard Model 
A+ s because, while the CKM suppression persists, we now have in addition a strong sup- 
pression from small Yukawa couplings. Consequently, there is no effect on Aqj> or on R. 

Neutral Higgs exchange in models without natural flavor conservation (NFC) could 
contribute to both b — > uus and b —>■ dds. But if the smallness of scalar mediated 
FCNC is explained by a horizontal symmetry, then we expect a suppression of order 
< 0(i7ibm U: d/m 2 z ) ~ 10~ 5 , which means that the effects are negligible. 

In left-right symmetric models there is a new contribution from W^-mediated decays, 
but it is suppressed by 0(m^ /mf^ ) < 10~ 2 . The CKM ratio is expected to be 0(1). So, 
again, we expect no observable effects on either Aqp or R. 

To summarize the situation regarding new physics effects, we note the following points: 

(i) There are several well-motivated extensions of the standard model which can signifi- 
cantly affect B — ► 71K decays. 

(ii) In models where there are new tree diagram contributions, there are no relations, in 
general, between the contributions to B + — > tt + K° and the contributions to B° — > ir~K + . 
The new physics effects in the amplitudes (|1.3|) introduce six new parameters. 

(Hi) In models where the new contributions are through QCD penguin or chromomag- 
netic dipole operators, isospin is a good symmetry, and the number of new parameters is 
therefore reduced from six to three. 



(iv) We note that if a large CP asymmetry ( |1.4|) is measured, that will invalidate the 
bound ( |1.9| ). In contrast, if a small CP asymmetry is measured, that would not provide 
an unambiguous confirmation for the validity of this bound, because it could be a result of 
small strong phases rather than a small magnitude of final state interactions. 
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V. CONCLUSIONS 

We have analyzed the effect of final state interactions on the search for New Physics in 
B — ► nK decays. Using a phenomenological model, we found that while such effects are 
unlikely to be large enough to dominate individual branching fractions, they can still compli- 
cate those avenues for identifying New Physics which rely on a Standard Model suppression 
of weak phases in the matrix elements for B ± — > k^K. As a result, and in contrast to pre- 
vious expectations, we conclude that the observation of Aqp ~ 0.2 or sin 2 7 ~ 1.2R would 
not be an unambiguous sign of a source of CP violation beyond the CKM matrix. While 
we do not claim to compute the magnitude of final state interactions reliably, our model is 
sufficient to demonstrate that effects of this size are entirely generic and cannot be ruled 
out without independent empirical evidence. We propose a simple test which could probe 
this question experimentally. Finally, in anticipation of the observation of CP violation in 
these channels at a level which cannot be explained by the Standard Model, we discuss the 
features of various models of New Physics as they would be manifested in these decays. 
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